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Thinning  of  a Point  Process  Over  Time 
by 

Richard  F.  Serfozo,  Syracuse  University 

1.  ^troduction 

The  first  study  of  thinnings  of  point  processes,  Renyi  (1956),  contains 

the  following  result.  Suppose  that  points  are  placed  on  the  nonnegative  real 

line  R_|_  such  that  the  interpoint  distances  are  independent  and  identically 

distributed  with  mean  c.  These  points  are  thinned  over  time  (or  in  stages) 

such  that  at  time  n each  point  that  has  survived  till  then  is  independently 

retained  with  probability  p^  and  deleted  with  probability  1 - p^.  Let 

denote  the  number  of  points  still  remaining  in  the  interval  [0,  where 

a = (p,,..p  ) ^ ->■  Think  of  N as  the  thinned  process  at  time  n with  as  a 
n *^1  *^n  n ' n 

new  unit  of  scale  for  R^.  Then  converges  in  distribution  as  n • to  a 
Poisson  process  with  intensity  c 

As  an  Illustration  of  this  result,  consider  a production  system  in  wliiili 
parts  (of  one  type)  are  produced  and  then  sent  througii  a series  of  n work 
stations.  Defective  parts  occur  according  to  a renewal  process,  .ind  I'.icli  i 
of  these  will  be  detected  and  discarded  witii  probabilitv  1 - p at  t lie  m-tli 

ID 

station.  Consequently,  the  flow  of  defective  parts  is  titinned  over  n stagi-s. 

If  n is  large  so  that  p^^...p^  is  very  small,  then  according  to  Renyi 's  nsiilt, 
the  flow  of  defective  parts  from  the  n-th  station  can  be  approximated  by  a 
Poisson  process. 

Since  1956  twenty-four  articles  (see  the  references)  have  been  written 
on  various  aspects  of  the  following. 


This  research  was  sponsored  in  part  by  the  Air  Force  Office  of  Scientific 
Research  under  Grant  AFOSR-74-2627 , and  by  the  National  Science  Foundation 
under  grant  ENG-75-13653. 


Thinning  Problem.  Points  are  placed  in  a space  according  to  a certain  probability 
law,  and  then  the  points  are  randomly  deleted  according  to  some  rule.  Is  it 
possible  to  normalize  the  thinned  process  by  a rescaling  of  the  space  (as  above 
with  the  a^)  so  that  it  converges  in  distribution  to  some  process  as  the  thinning 
becomes  complete?  And  what  are  the  possible  norming  constants  and  limits? 

This  problem  is  analogous  to  the  central  limit  problem  for  sums  of  random 
variables,  which  is  concerned  with  finding  norming  constants  and  such  that 
^ ^k~l^k~^n  related  stochastic  process  [3,  Theorem  10.1])  converges  in 

distribution.  Here  the  or  its  variance  is  converging  to  infinity,  and 

in  thinning,  the  thinned  process  is  converging  to  zero.  The  rescaling  by  B^ 
is  comparable  to  rescaling  the  space  in  a thinning.  The  central  limit  theory 
has  yielded  many  insights  into  random  phenomena  in  nature  that  arise  as  sums  of 
a large  number  of  random  variables.  For  example,  we  can  now  readily  identify 
various  types  of  random  phenomena  that  can  be  modeled  as  (or  approximated  by) 
normal  or  infinitely  divisible  random  variables,  Weiner  processes,  stable 
processes,  etc.  A knowledge  of  thinning  will  yield  similar  insights  into  point 
processes  that  arise  from  thinning  procedures  such  as  the  following:  (1)  debugging 

of  complex  computer  programs,  weapons  systems,  communications  networks  etc. 

(a  check-out  routine  describes  a random  path  through  the  system  and  the  residual 
errors  form  the  thinned  process),  (2)  filtration  of  pollutants  in  the  air  or  water 
by  magnetized  smoke  stacks,  car  mufflers  or  charcoal  beds,  (3)  elimination  of 
undesirable  insect  or  animal  pests,  plants  or  cell  growth,  and  (4)  cutting  down  on 
potential  disease  carriers  by  an  immunization  program. 

In  a recent  study  [28]  we  showed  how  many  thinnings  of  point  processes  can 
be  analyzed  in  terms  of  compositions  and  inverses  of  random  measures.  In  doing 
so  we  unified  many  of  the  results  on  thinnings.  The  main  emphasis  in  [28],  similar 
to  most  of  the  studies  referenced  herein,  was  on  a single,  very  thorough,  thinning 


operation.  In  this  article  we  consider  thinnings  in  which  points  are  deleted 
over  time  or  by  a sequence  of  thinning  operations.  In  Section  3 we  show  that 
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the  asymptotic  behavior  of  a thinned  process  is  essentially  equivalent  to  that 
of  the  extreme  values  (high  level  exceedances)  of  the  lives  of  its  points  under 
the  thinning.  We  use  this  result  in  Section  4 to  describe  independent,  regenerative, 
and  semi-stationary  thinnings  of  point  processes  on  R_^.  Then  in  Section  5 we 
describe  independent  thinnings  on  more  general  spaces. 

2.  Notation. 

We  shall  use  the  following  notation  for  point  processes  and  random  measures. 

This  is  similar  to  that  in  the  basic  works  [1],  [9]  and  [13]. 

Let  M denote  the  set  of  measures  on  theBorel  sets  B_|_  of  R^.  that  are  finite 

on  bounded  sets.  We  endow  M with  the  vague  topology  [1].  This  topology  is 

metrlzable  and  the  following  are  equivalent  statements. 

(i)  the  vague  topology. 

(ii)  /f  (x)  dy^(x)->-/f  (x)d|j  (x)  for  each  continuous  f on  R^  with  compact  support, 

(iii)  u^(a,b]-»-u(a,b]  for  all  a,  b such  that  y({a,b})  = 0. 

Let  M be  the  smallest  o-field  on  M containing  the  open  sets  of  the  vague  topology. 
This  is  the  same  as  the  smallest  o-field  tliat  makes  the  mappings  y->-u(A),  for 
Ar 8^,  measurable. 

A random  measure  K on  R_j_  is  defined  to  be  a measurable  mapping  from  a 
probability  space  to  (M,M) . We  denote  by  C(A)  the  random  variable  describing  the 
mass  in  the  set  AzB^,  and  we  let  C(t)  = 5(f0,t])  for  ttR_|^.  If  L(t)  is  integer- 
valued for  each  t,  then  C is  called  a point  process. 


A sequence  of  random  measures  converges  in  distribution  to  a random 

measure  C,  written  1 5r  if  the  distribution  of  I converges  weaklv  to  the  distribution 

of  £.  That  is,  Eh(C  ) ->Eh(C)  for  each  bounded  continuous  function  h on  M,  sec  [3]. 

n 
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Tfie  following  are  equivalent  statements: 

(1)  c„Pt. 

(ii)  j f(x)d^^(x)  / f(x)d^(x)  for  each  continuous  f on  with  compact  support. 

V 

(iii)  , . . . , C^(Aj^) ) > (5(A^),...,  C(A|^))  for  any  bounded  intervals  A^,..., 

Aj^  in  R_^  satisfying  C(<SAj^)  = ...  = C(6Aj^)  = 0 a.s.,  where  6A  denotes  the  boundary 
of  A. 

The  above  terminology  also  holds  with  obvious  modifications  for  measures  on 
locally  compact  second  countable  Hausdorff  spaces  such  as  R*^. 

Much  of  our  analysis  will  deal  with  compositions  of  measures.  If  ^ and  n 
are  random  measures  on  R_^  then  their  composition  ^oq  is  a random  measure  witli 
^on(t)  = ^(n(t)).  We  also  let  ^oc,  for  ceR_|_,  denote  the  measure  with  ^oc(t)  = ^(ct)  , 
In  the  following  r-  suit  the  C's  and  q's  are  random  measures  on  R^  and  the  a's 
and  B's  are  random  variables. 

Lemma  2.1.  (i)  If  ,r\  ) ^ (Cjd).  where  E(t)  is  continuous  a.s.  or  n(t)  is 

n n 

strictly  increasing  and  n(0)  = 0 a.s.,  then  ^ Cop. 

(ii)  If  (C  ,a  >B  ) ^ where  a <6  and  ^({a,6}  ) = 0 a.s.,  then 

n n n n n 

Proof . These  assertions  follow  from  Theorem  3.2  and  Corollary  2.4  in  [28]. 

3.  Thinning  High-Leyel  Exceedances  of  Stochastic  Processes 

The  thinning  procedure  that  we  shall  consider  here  and  in  the  next  section 
is  as  follows.  Points  are  placed  in  R^  at  random  sites  0 5 < S,  < ...  whcrt* 

Sj^  >■  “>  a.s.  Throughout  Sections  3 and  4 we  assume  the  following. 

-1  V 

General  Assumption,  k Sj^  ‘ c where  c is  a positive  constant. 

The  points  are  deleted  over  time  (or  in  stages)  according  to  some  rule.  We 


Jet  Yj^  denote  the  random  length  of  time  - which  we  take  to  be  discrete  - that  tlic 
point  at  site  Sj^  survives  under  the  thinning.  (In  this  section  we  make  no 
assumptions  on  the  dependency  between  the  Y's  or  S's.)  Then  at  time  n,  the 
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number  of  points  still  remaining  in  a set  At8^  is  given  by 


= I I(Yj^  > n)6g  (A) 
k k 


where  I(Yj^  > n)  = 1 or  0 according  as  Y^^  > n or  n,  and  ^g(‘)  denotes  the  Dirac 

measure  with  unit  mass  at  s.  In  this  section  we  show  how  the  asymptotic  behavior 
of  the  thinned  process  is  related  to  large  values  of  the  Y's.  In  the  next 

section  we  present  more  specific  results  when  the  Y's  are  a regenerative  or  semi- 
stationary sequence. 


To  record  large  values  of  the  Y's  we  shall  use  the  process 
[t] 


= I I(Y  > n) 
k=l 


where  [t]  is  the  integer  part  of  t.  The  is  the  exceedance  process  at  level  n 
associated  with  the  Y's.  Namely,  is  the  number  of  the  lifetimes  Y^,...,Y^^ 

which  exceed  n. 

The  following  result  shows  that  the  convergence  of  is  equivalent  to  the 
convergence  of  n^.  Note  that  the  randomness  of  the  S's  affects  the  convergence 
of  only  through  the  constant  c.  Here  and  below  the  a^  and  b^  are  positive 
constants  with  a <»,  and  f and  n are  random  measures. 

Theorem  3.1.  A necessary  and  sufficient  condition  for  b oa  some  C is  that 

-IV  V 

b n oa  some  n.  In  either  case  C = noc. 
n n n 

Proof . Clearly 

C(t) 

C^(t)  I(Y^  > n)  = n^(C(t)) 

k=l 

where  /;(t)  = max  {k  : Sj^<  tl.  Also  note  that  the  assumption  k c is 

-1  P -1 

equivalent  to  t c(t)  >■  c . Thus  the  assertion  follows  by  [28,  Theorem  4.1]. 


Rema r k 3.2.  Theorem  3.1  also  holds  If  convergence  in  V is  replaced  throughout  by 
convergence  a.s.,  or  convergence  in  probability.  This  follows  from  ttie  nature  of 
its  proof. 


4.  S^tiona^,  Regenerative,  and  Semi-stationary  ThJ.nnings 

In  this  section  we  describe  the  asymptotic  behavior  of  the  thinned  process 
5 =5^  > n)6„  , as  described  above,  for  several  types  of  thinnings.  We 

label  these  thinnings  according  to  the  dependency  structure  on  the  Vs. 

We  begin  with  an  independent  thinning. 

Theorem  4.1.  Suppose  Y2,  ...  are  Independent  with  a common  distribution.  Then 

V -1 

C^oa^  ->•  C,  where  C is  a Poisson  process  with  intensity  Ac  , if  and  only  if 
anP(Yi  > n)  ->  A > 0. 

Proof . This  is  a special  case  of  Theorem  4.4  below. 

Examgl^4.2.  (Renyi's  result)  Suppose  the  thinning  is  such  that  at  time  n 
each  point  that  has  survived  is  independently  deleted  with  probability  1 - p^ . 

That  is  P(Y,  > n)  = p,  ...  p . If  a = (p  . . . p ) ^ tlien  1 oa  converges  in 

distribution  to  a Poisson  process  with  intensity  c 

Our  next  result  is  for  stationary  thinnings.  In  it  we  use  the  following. 
I)e^ndency_^  C^pndl^ions  4.3.  Let  (k,  n > 1)  be  nonnegative  integer-valued  random 


variables . 


(i)  For  any  integers  1 < < k2  •••  < k^  < ...  < > with 

I = •••  = \k  = ^n£,  = = ^ne  = 0^ 

1 pi  q 


-P(X 


nk. 


= X , = 0)  P(X  , = . . . = X „ = 0)  d a . 

nk  nl,  ni  ' ni 

p ] q 


where  a . is  nonincreasing  in  1 and  a . -*-0  for  some  sequence  1 -»■ 

ni  ni  n 

n 


- k > i 
1 P 


<”  such  that 
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(li)  lim  sup  a ) P(X  , > 1,  X , > 1)  = o(m  as  m 
n n > mn , 1 ~ ’ mn , j 

j=2 


These  dependency  conditions  are  essentially  those  used  in  [16]  and  [17] 

(the  n and  u^  in  these  references  are  our  a^  and  n respectively).  Note  that 

( i)  is  slightly  weaker  than  the  usual  mixing  conditions  that  describe  dependenc it s . 

Suppose  that  Y^,  Y2,  ...  is  a strictly  stationary  sequence  such  that  tin. 

X^j^  = I(Yj^  > n)  satisfy  Conditions  4.3  and 

a P(Y,  > n)  ->■  A > 0. 
n 1 

V -1 

Then  C where  ^ is  a Poisson  process  with  intensity  Ac 

Proof . This  is  a special  case  of  Theorem  4.7  below. 

We  now  consider  regenerative  thinnings.  Namely,  we  shall  assume  that 
Yj^,  Y2,...  is  a regenerative  process  over  the  integer-valued  indices 
1 ~ "‘’g  ^ T2...  That  is,  the  vectors 

<'k  - ’k-r  “•  “••••>  f"rkki 

are  Independent  and  identically  distributed.  This  might  occur  when  the  points 
are  thinned  in  independent  groups  of  random  sizes.  The  above  vector  then  would 
describe  the  lives  of  the  points  in  the  k-Lh  group  which  consists  of  t 
points. 

As  an  example,  suppose  the  initial  points  at  sites  S^,  S^,  ...  have  generic 
makeups  whicli  determine  their  lifetimes  Yj^,  Y2,  ...  such  that 

P(Y|^  = n I = i.  Y^(?,  k;)=  P(Y^  = n | = i)  . 

Assume  that  is  a positive  recurrent  Markov  chain  with  a countable  state  space 
and  = i.  Let  1 = *^0  ^ successive  indices  k for  which  B]^  = >• 

Then  Yj^  is  a regenerative  process  over  the  I's.  Another  interpretation  of 
this  example  is  that  the  initial  points  arc  all  the  same,  but  the  thinning  is 
inhomogeneous  in  that  at  site  Sj^  it  depends  on  tlie  parameter  Bj^. 
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For  the  next  result  we  assume  tliat  Y,^, 


is  a regenerative  process 


over  I = Tq  < < ...  with  u = K(ij^  - Tq)  < and  we  let 


il-l 

I ’ n) 

k=l 


for  k > I. 


We  also  let  f;  be  a random  measure  such  that  f,(t)  has  stationary  independent 
increments  with 

-sf.(l).  ^ 


E(e 


^OO 

) = exp-  f (1 


e ®^)dp(x) , 


where  [j  is  a measure  on  (0,  “)  satisfying  min  {1,  x}dp(x)  < «>. 

Theorem  4.5,  For  b oa  > fo(ac)  ^ it  is  necessary  and  sufficient  that  the 
following  hold. 

(1)  a P(b  < t)  > p(t)  for  all  t with  ;if{  t r)  = 0. 

n n n 


(2)  lim  lim  a b 

6.->U  n n n 

n K» 


lim  lim  a b ^ EX^ 
e ^0  n ■■■”  ti  n n 


0,  where  X is  X truncated  at  l. 
n n 


-1, 


Remark  4.6.  Note  that  b^  ^n'^^n  distribution  to  a compound  Poisson 

process  whose  jumps  occur  according  to  a Poisson  process  with  intensity  (tc)  ^ and 
whose  jump  sizes  have  distribution  u if  and  only  if  (1)  holds  with  u(0.°-)  = 1. 

(The  latter  implies  (2).) 

Pri!of.  Let 


-1 


r.  (t)  = b ' V"  X 
" " k=l 


, , where  X , 
nk  nk 


I (Y 


n)  . 


k-1 


Conditions  (1)  and  (2)  are  necessary  and  sufficient  for  t (1)  ■ f.  ( 1 ) , see 

[13,  Theorem  6.1]  which  is  a generalization  of  [4,  p.  564].  Also  '^^(1)  * '(1) 

is  equivalent  to  t ? f,  see  [5,  p.  480].  In  adaition,  h ^ t,  oa  ^ fo(ac)  ' is 
^ n n n n , , 

-1  P -1  V 

equivalent  to  b n oa  ^ ''.on  by  Theorem  3.1.  Recall  that  n (t)  = ' 

' n n n n 

Thus  to  prove  the  theorem, it  remains  to  show  that 

(4.1)  b n oa  > \o'i  <=>  4 • • 

n n n ^ ^ ^n 


k=l 


n } . 


To  this  end  we  shall  use  the  representation 
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(4.2)  b oa  = 4 oy  + (b  oa  --  4,  OY„)  . 

nnn  n n nnn  nn 

whore 

y (t)  = a ^v(a  t)  and  v(t)  = max  (k  : x,  £ tl. 
nnn  ^ 

Wo  shall  also  use  the  fact  that 

(4.3)  ■>  a~^A  where  A is  the  Lebesque  measure. 

This  follows  since  “ "'^1^  1 independent  with  a common  distribution  and 

j -1  /.X  -1 

mean  «,  and  t v(t)->-a  a.s. 

To  prove  (4.1)  we  first  assume  that  4 ^ 4.  From  this  and  (4.3),  it 

n 


follows  by  Lemma  2.1  that 

//  /X  ^ r -1a  a-  -1 

(4.4)  4 oy  -*  4oa  A = f.oa  , 

n n 


and 

(4.5) 


-1  -K  V 

b n oa  - 4 oy  <4  o(y  + a ) - 4 oy  -*-0. 
n nn  n n-n  n n n n 


(Here  4 o(y  + a ^)(t)  = 4 (y  (t)  + a ^).)  Then  from  (4.2)  - (4.5)  it  follows  that 
nnn  n n n 

K-1  i’  r -1 

b n oa  -F  Aoa 
nnn 

Now  assume  the  latter  holds.  Clearly 


4 = 4 oy  oy  = (b  n oa  0({)  )oy 

n nnn  nnn  n n 

-1 


whore  i (t)  = a r . . and 

n n v(a  t) 

n 


Y^(t)  = inf  {u  : y^(u)  > t)  = r 

n 


The  latter  follows  by  standard  properties  of  inverses.  Under  our  assumption  we  have 
(b^^n^oa^,  (+>^,  y^)  ^ (40a  \ A,  ah). 


Tlu-n  bv  Lemma  2. 

V 


-1 


4 ► foa  oaA  = 4- 

n 


Til  is  eompletos  the  proof. 
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We  now  consider  semi-stationary  thinnings.  Namely,  we  assume  that  is  a 
semi-stationary  process  over  the  integer-valued  random  variables  1 ~ ^0  ^ ^2’’ 

This  means  that  the  sequence  of  vectors 

_i.0,0,...)  for  k > 1 

k-1  k k 


' I 


is  stricfly  stationary,  see  [27].  Special  cases  of  these  processes  are  stationary 
processes,  positive  recurrent  Markov  chains,  regenerative  processes  and  many 
combinations  of  these  types  of  processes. 

As  an  example  suppose  that  the  point  at  site  has  a generic  makeup  described 
by  and  the  thinning  operation  at  is  based  on  a parameter  such  that 

\ 

where  f is  some  function.  Suppose  is  a positive  recurrent  Markov  chain  on  a 

countable  state  space  and  1 = < . . . are  the  successive  indices  k for  which 

6j^  = i.  And  suppose  that  is  a stationary  sequence  which  is  independent  of 
Then  is  semi-stationary  over  the  x's. 

The  next  result  on  semi-stationary  thinnings  contains  Theorems  4.1  and  4.4 
as  special  cases. 


Theorern  4.7.  Suppose  Yj^,Y.,,...  is  semi-stationary  over  1 = x^  < x^  < ...  with 

I 

"=Vi 


ct  =(Ex,  - Tr,)  ^ " such  that  X , 
i U nk 


l(Yj^  > n)  satisfies  Conditions  4.3  and 


^ j, 


a P(X  = 1)  7 > 0,  and  a P(X  , i 2)  0. 

n nl  n nl 

Then  t oa  ^ E where  C is  a Poisson  process  with  intensity  l(uc) 
Ti  n 


Proof.  Let 


where 


■ J,  \k  - 

k=  I 
[a^t] 

V')  - I 

k=  1 


u 


From  [17,  Theorem  3.2]  It  follows  that  Q.  Also  for  each  t 


[at] 

P(^„(t)  - C(t)  = 0)  = P(  n < It) 


n 

[a„t 


k=l 


= ' - ^ 2})  > 1 - [a„t]P(X^^  > 2)  ^ 1. 


That  is,  ^ 0.  We  then  have  ^ C.  From  this  it  follows,  by  an  argument 

as  in  the  second  and  third  paragraphs  in  the  last  proof,  that  Hj^oa^  Coc.  But 


this  is  equivalent  to  C^oa^  ->■  C by  Theorem  3.1. 


5.  Indeg^^n^Thinning  in  General  Spaces. 

In  this  section  we  present  two  analogues  to  Theorem  4.1  for  independent 
thinning  in  a locally  compact  second  countable  Hausdorff  space  X such  as  R*'. 

The  Theorems  4.4,  4.5  and  4.7  do  not  have  such  clear-cut  analogues,  since  the 
notions  of  regeneration  and  semi-stationar ity  require  a total  ordering  of  the 
initial  points  which  can  be  specified  in  many  ways. 

For  our  first  result,  we  assume  that  points  are  placed  in  X at  random  sites 
Sj^,S2,...  and  are  deleted  over  time  such  that  their  lifetimes  Y^,Y2,...  are 
independent  and  identically  distributed  and  are  independent  of  the  S's.  At 
time  n the  thinned  process  is 

■ 1 

k 

This  C on  X has  the  same  structure  as  before.  We  also  let 
n 

? = f <5  -Ip  and  f.  oa  (A)  = C [a  x : xeA}. 

n R a S,  n n n n 

n k 

00  , 


where  a 


]2 


Theorem  5.1.  If  a ?y,  where  y is  a measure  on  X,  and  a„P(Y,  > n)  -/  X > 0, 
— ~~ — n n n i 

p 

tnen  £ oa  > E,  where  £ is  a Poisson  process  with  intensity  measure  Xy. 
n n 

Proof.  For  any  bounded  disjoint  Borel  sets  A , . . . ,A,  in  ^ have. 


^n°^n^V^  = < ‘ 

where  . . . . , are  independent  copies  of  ^ also 

independent  of  i;^.  Under  the  assumptions 

(n^lOa^,...,nnk°V  ^ 

where  are  independent  Poisson  processes  with  intensity  A (a  constant). 

We  also  have 

^^n^^n(Ai).---.a;\„(Aj^))  - (y(Ap....,y(A^)). 

Then  it  follows  by  Lemma  2.1  that 

(£n°^n^^l^  ’ ■ ■ ■ ’ ^n°^n^^k^^  ^ ( ni(y  (Aj^)  ),...,  rikCh  (Aj^)  ) = (£ (A^^)  ,...,£  (Aj^)  ) . 
This  completes  the  proof. 

The  next  result  is  a slight  variation  of  the  proceeding  and  its  proof  is 
the  same.  For  this  we  assume  that 

«„  - I =■  ">‘s  ^ ■ ”s  , 

nk  k nk 

where  S , is  the  site  of  the  k-th  point  at  time  n and  the  Y's  are  as  above.  This 
nk 

describes  thinnings  in  which  the  points  may  move  over  time.  We  also  let 
5 be  a random  measure  on  X. 


Theorem  5.2.  If  a^Z,  ^ t,  and  a P(Y,  > n)  •+  X > 0,  then  £ ^ £,  where  £ is  a 

— n n n 1 n 

conditional  Poisson  process  with  random  intensity  measure  X£. 


H'.  ^ I yif'^nraMrri  r i'  'i 


The  conditional  (or  doubly  stochastic)  Poisson  process  C with  intensity 
\r,  is  defined  by 


P(C(Aj)  = n^,  . . . ,C(A|.)  = n^) 

= e ^..e  '"/(n^!  . . .n^!)P(A^(Aj)edu,  . . . , Xf.(A^)£dUj^)  . 

These  processes  were  first  studied  by  Cox  and  Lewis,  see  [19]  and  [26]  and  theii 
references . 

Theorem  5.2  for  the  Renyi  thinning,  where  P(Y^  > n)  = p’  , was  proved  in 
[ 19 ] and  in  [ 12] . 
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